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ABSTRACT 

We investigate whether gamma-ray pulsars viewed at a large angle to the neutron star magnetic pole 
could contribute to the new population of galactic unidentified EGRET sources associated with the 
Gould Belt. The faint, soft nature of these sources is distinctly different from both the properties of 
unidentified EGRET sources along the galactic plane and of the known gamma-ray pulsars. We explore 
the possibility, within the polar cap model, that some of these sources are emission from pulsars seen 
at lines of sight that miss both the bright gamma-ray cone beams and the radio beam. The off-beam 
gamma-rays come from high-altitude curvature emission of primary particles, are radiated over a large 
solid angle and have a much softer spectrum than that of the main beams. We estimate that the 
detectability of such off-beam emission is about a factor of 4-5 higher than that of the on-beam emission. 
At least some of the radio-quiet Gould Belt sources detected by EGRET could therefore be such off-beam 
gamma-ray pulsars. GLAST should be able to detect pulsations in most of these sources. 
Subject headings: stars: neutron; pulsar: general; gamma rays: pulsars 



1. INTRODUCTION 

The origin of the 170 unidentified gamma-ray point 
sources detected by the Energetic Gamma-Ray Experi- 
ment Telescope (EGRET) on board the Compton Gamma- 
Ray Observatory (CGRO) throughout its mission (Hart- 
man et al. 1999) is one of the most intriguing questions of 
high-energy astrophysics. Grenier & Perrot (1999) found 
that some of the unidentified sources are significantly cor- 
related with the Gould Belt of massive stars, a nearby 
galactic structure surrounding the Sun. The Gould Belt is 
an expanding disk of gas and young stars, most with ages 
less than 30 million years, inclined about 20° to the galac- 
tic plane. Recently, Gehrels et al. (2000) studied the spa- 
tial and flux distribution of the 120 steady sources in the 
third EGRET catalog and found that they divide into two 
groups: higher flux sources distributed along the galactic 
plane with latitudes less than 5°, and a new population of 
lower flux sources at mid-latitude which indeed seem to be 
associated with the Gould Belt, at a relatively nearby dis- 
tance of 100-300 pc. The two groups also show significant 
differences in their spectra, with the low-latitude sources 
having average photon spectral index of 2.18 ± 0.04 and 
the mid-latitude sources having average photon index of 
2.49 ± 0.04. None of these sources have known counter- 
parts at other wavelengths. Grenier & Perrot (1999) sug- 
gested that the Gould Belt 7-ray sources could be young 
pulsars, formed in supernova explosions of massive stars 
of the belt. However, the relatively soft spectra and lu- 
minosities of these sources, averaging ~ 6 x 10 30 ergs 



are not characteristic of the known 7-ray pulsars detected 
by EGRET, which have spectral indices of 1.5 — 2.0, lumi- 



the bright, hard 7-ray beams but detecting only the off- 
beam emission. According to the polar cap model (e.g. 
Daughcrty & Harding 1996 [DH96]), 7-ray emission occurs 
throughout the entire pulse phase. A synchrotron-pair cas- 
cade at low altitude radiates the hard on-beam emission 
in a hollow cone centered on the magnetic pole to pro- 
duce the bright double-peaked pulses. Primary electrons 
continue to radiate curvature emission on open field lines 
to high altitudes beyond the cascade region, producing a 
lower level of softer off-beam emission. Due to the flaring 
of the dipole field lines, this emission may be seen over 
a large solid angle, far exceeding that of the main beams. 
Since the radio emission is expected to originate within ten 
stellar radii of the neutron star surface, it is quite prob- 
able to see off-beam 7-ray emission and miss the radio 
beam. Such off-beam or off-pulse emission may have been 
detected by EGRET in the Crab, Vela and Geminga 7-ray 
pulsars during non-pulse phases (Fierro et al. 1998). We 
show that this expected off-beam emission has the right 
characteristics (spectral index, luminosity and numbers) 
to account for the EGRET sources associated with the 
Gould Belt. 

2. OFF-BEAM PULSAR EMISSION 

2.1. Geometry 

The basic geometry of polar cap emission, projected 
onto the stellar surface, is shown in Figure 1. The main 
beam of the 7-ray cascade emission zone is located along 
the edge of the hollow cone with opening angle 7 , cen- 



tered on the magnetic axis, denoted by the dotted circle 
labeled 7. The altitude of this zone is about 1-2 stellar 
radii above the stellar surface (DH96, Harding & Mus- 
limov 1998). The off-beam 7-ray emission zone is located 
at altitudes above the cascade emission zone, from several 
stellar radii to the light cylinder. Off-beam emission visi- 
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nositics 10 — 10 ergs and 0.5-4 kpc distances. 

Here, we show that the bulk of the Gould-Belt 7-ray 
sources could be pulsars if their emission is seen at large 
angles to their magnetic axes, such that we are missing 
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ble to EGRET above 20 MeV is shown as a shaded region. 
The radio emission can have a core component, thought to 
be located at low altitude along the magnetic axis, and/or 
a conal component (s), thought to be located in a hollow 
cone at altitudes of tens of stellar radii (e.g. Rankin 1993). 
In the polar cap model, the radio emission originates from 
the secondary pair cascade. In Fig. 1, the possible lo- 
cation of radio conal emission with opening angle 9r is 
shown as a dashed circle labeled R. The circles centered 
on the rotation axis, labeled £i and (2, denote the sweep 
in pulse phase of different observer lines-of-sight viewing 
on-beam and off-beam 7-ray emission respectively. 



b) 




Fig. 1. — The geometry of the polar cap model projected onto the 
stellar surface. The dotted circle labeled 7 is the main gamma-ray 
conal beam, the dashed circle is a conal component labeled R of 
the radio emission beam, and the circle enclosing the shaded region 
denotes the cone where the typical curvature radiation emission of 
the primary particles is 100 MeV. The shaded area is the above-100 
MeV-emission area of the pulsar. The inclination angle between f2 
and B is a and £'s denote different lines of sight. The line-of-sight 
£2 cuts off the main gamma-ray beam and misses the radio emission 
beam, and we identify such a configurations as the off-beam pulsar 
candidates for the "Gould Belt" faint, soft, unidentified EGRET 
sources, (a) a not small: £1 gives a typical Vela-type emission con- 
figuration; (b) a small: £1 gives a typical Gcminga-type emission 
configuration. 



The off-beam gamma-ray emission is mainly produced 
by curvature radiation of the primary particles, well be- 
yond the acceleration and cascade zone. For a certain 
inclination angle a and line of sight £, the off-beam emis- 
sion zone, between roB(niin) and roB(max), (or the cor- 
responding #oB(min) and 6*oB(max)) may be determined 
by geometry and the curvature radiation properties. The 
minimum height, roB(min), is the lowest height visible to 
a given line-of-sight and is defined by the radius at which 
the tangent to the last open field line, which is |6> for a 
dipole, where 9 is the polar angle, points in the direction of 
the observer (which is also the outer edge of the off-beam 
7-ray emission cone at that radius): 



r OB (min) ~ - 



!«-«) 



(i) 



For P — 0.3 s, and £ — a = 10°, we have r OB (min) ~ 19i?. 

The maximum height at which EGRET-band emission 
(> 100 MeV) can be produced may be defined as the ra- 
dius roB( max ) along the last open field line where the 
curvature radiation photon energy Eq = 100 MeV. The 
typical curvature photon energy is Eq = (3/2)(7 3 ?ic/ p) = 
322MeV7fP" 1 / 2 (r/i?)- 1 / 2 for dipolar field configuration, 
where 77 = 7/IO 7 is the primary particle Lorentz factor, 
p is the magnetic field radius of curvature, P is the pulsar 
period and r and R are the radii of the emission point and 
stellar radius, respectively. To be detected by EGRET, 
one thus requires 

77 > 77(100 MeV) = 7100 = 0.68P 1/6 (r/i?) 1/6 . (2) 

The expression for the Lorentz factor of a particle trav- 
eling along an open dipole field line and subject only to 
curvature radiation losses is (Harding 1981) 



To 



8 r R 



-1/3 



(3) 



where r e = e 2 /mc 2 is the classical electron radius and 70 
is the initial Lorentz factor at Rq, the radius of the accel- 
eration zone. Using Eqn (|3|), Eqn (0), and the expression 
for the last open dipole field line, 9 1 = (rfl/c), we have 
the following equation for roB(niax), 



7ioo (r OB (max)) 



7o,7 



g!in 7 gM r °» (max) ) (4) 

8c 1 R 1 w 



where 70,7 = 70/IO 7 . Finding the root of this equation 
yields values of roB(max) for different values of SI = 2ir/P, 
70 and Rq, which will be limited by the light cylinder dis- 
tance, ~ tlc I sin a in short-period pulsars. For P = 0.3 
s, 70 = 2 x 10 7 and R = 2R, one finds roB(max) = A8R. 

2.2. Relative Observability 

The relative opening angle, 9 OB ~ (fir OB (max)/c) 1 / 2 , of 
the "100 MeV cone" with respect to the on-beam 7-ray 
cone, 7 ~ (fh^/c) 1 / 2 , is: 9 OB /9 7 ~ (r OB (max)/r 7 ) 1 ^ 2 . 
For a fixed a, the relative detectability of the off-beam 
pulsars and the gamma-radio plus Geminga-like pulsars is 
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(note that only geometric effects, but no luminosity selec- 
tion effect is taken into account): 



cos(Max[0, a ■ 



,]) - cos(Min[7r/2, a + 6 OB ]) 



~ cos(Max[0, a - 7 ]) - cos(Min[7r/2, a + 6» 7 ]) 

Assuming a random distribution of a, the relative de- 
tectability of the off-beam pulsars and the gamma-radio 
plus Geminga-like pulsars is (Emmering & Chevalier 
1989): 

D = (1 - cosfl OB ) + (tt/2 - fl OB )sinfl OB _ 

(l-cos0 7 ) + (7r/2-0 7 )sin0 7 ' [1 

Given typical values P = 0.3 s, r 7 = 2, and assuming emis- 
sion along the "last open field line" , one gets OB /0 7 ~ 5 
and D ~ 4 according to eq. . Figure 2 shows the depen- 
dence of both r OB (max) and P/ on pulsar period. If we 
regard Geminga and PSR 0656+14 (the strongest gamma 
ray pulsar candidate) as the on-beam gamma-ray pulsars 
in the Gould Belt, we expect that there might about 8 
off-beam gamma-ray pulsars, a sizable portion of the 20- 
40 Gould Belt unidentified EGRET sources (Grenier et 
al. 2000). The luminosity selection effects will modify 
this number, but the above estimate suggests that at least 
some of the Gould Belt sources could be off-beam 7-ray 
pulsars. 
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Fig. 2. — Solutions for t 
sion, r OB (max), from Eqn. ( 

relative to on-beam emission, D, Eqn (&), as a function of period 



Figure la and lb show two different possibilities for ob- 
server orientation relative to the spin and magnetic axes. 
In Fig. la, £1 > a so that the line-of-sight cuts outside 
magnetic axes, seeing both bright rims of the 7-ray cone 
(and thus a double-peaked 7-ray profile) and the radio 
cone. Assuming the trailing radio component is missing, 
the broad band pulse profiles of the Vela pulsar and PSR 
B1046-58, with 7-ray pulse phase separation W = 0.4 
and phase difference between the radio and leading 7-ray 
pulses of 5<f) — 0.12, are consistent with this picture over 
a large a — Q phase space as long as a is not too large. 
In Fig. lb, Ci < a so that the observer cuts between 
the spin and magnetic axes, seeing both bright rims of 



the 7-ray cone (and thus a double-peaked 7-ray profile), 
but missing or grazing the radio cone. Such a configu- 
ration might describe Geminga-likc pulsars. Thus, the 
phase space Max(0, a — 7 ) < £ < Min(7r/2, a + 7 ) com- 
prises radio-loud 7-ray pulsars and Geminga-like pulsars, 
both with double-peaked 7-ray profiles. The phase space 
a + 8 7 <C<7r/2orO<C<Qi — 7 contains the off-beam 
7-ray pulsars having single peaked 7-ray profiles. Thus, 
we expect that the majority of the 7-ray pulsars detected 
in the Gould Belt will be radio quiet. 

2.3. Luminosity and Spectrum 
Since the energy loss rate for curvature radiation, 7, 

7 = f-V = 6.6 x l0^P~\r/Rr\ (7) 



declines with height, the luminosities of the off-beam 
sources then mainly depend on 7[r OB (min)]. A rough es- 
timate of the relative luminosities of the on-beam and off- 
beam 7-ray emission is then 



L 7 (on) 7(on) 
L 7 (off) ~ 7(off) 



To 



i -1 



7[r OB (min)] 



r OB (min) 



(8) 



The maximum Lorentz factor of the primaries is (Zhang 
& Harding 2000) 



7o 



1.4 x lO'p-^cosa) 1 / 4 , 



which is 70 7 ~ 1.9 for typical values. Using Eqn 
d) and Eqn®, for P = 0.3, 



L 7 (on) 
L 7 (off) 



80, 



(9) 
Eqn 

(10) 



consistent with the typical luminosity of the Gould Belt 
sources. 

Polar cap cascade simulations generally confirm the ap- 
proximate results derived above. A simulation run for the 
parameters of the Vela pulsar, using the cascade code de- 
veloped by DH96, give the brightness and model spectra 
of 7-ray emission as a function of £ — a, shown in Fig- 
ures 3 and 4 respectively. The contrast between on-beam 
and off-beam emission will be smaller for longer period 
pulsars. Small values of C — ol cut through the bright syn- 
chrotron core of the cascade, where the spectrum is a hard 
power law of index 2 with a high-energy cutoff above 
1 GeV due to electron-positron pair production. For in- 
creasing values of £ — a the line-of-sight falls off the main 
beam, views higher-altitude emission and the luminosity 
decreases. However, the tail of high-altitude cascade cur- 
vature emission extends to large values of £ — a and allows 
the off-beam pulsars to be visible over a large solid an- 
gle. Figure 4 shows that the off-beam emission spectrum is 
much softer in the EGRET range, because the high-energy 
cutoff falls well below 1 GeV. The energy of the cutoff de- 
creases with increasing £ — a, as r OB (min) increases and 
critical curvature photon energy Eq , which now defines the 
cutoff, decreases. The spectrum of the off-beam emission 
in the EGRET band appears very soft because the Eq is 
around 100 MeV. Power law fits to the off-beam model 
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spectra between 0.1 and 1 GeV (an example is shown in 
Figure 3) give photon indices in the range 2.2-3.0, consis- 
tent with what is observed for the Gould Belt sources. 




10 20 30 40 50 
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Fig. 3. — Distribution of 7-ray emission at 100 MeV and 1 GeV as 
a function of polar angle to the magnetic axis, f — a for parameters 
of a Vela- type pulsar (cf. Fig. 4). 

Given the relative luminosity of off-beam pulsars calcu- 
lated from Eqns (Q- |l^) , we can estimate absolute luminos- 
ity and the number of sources detectable by EGRET and 
GLAST. Using the predicted 7-ray luminosity, £ 7 (on) ~ 
5xl0 32 ergs -1 , of a pulsar with P = 0.3 s and B = 2x 10 12 
G from a polar cap cascade (Zhang & Harding 2000, 
Eqn [60]), an average luminosity for an off-beam pul- 



6 x 10 30 ergs 1 



With an out- 



sar would be I/ 7 (off) 
of-plane limiting point source sensitivity of § EGRET = 
6 x 10~ 8 phcm~ 2 s _1 (Hartman et al. 1999), we estimate 
that EGRET can detect off-beam pulsars out to a limiting 
distance of dy lm ~ 260 pc(fioB) -1 ^ 2 - Using Eqn (g) and 
the ratio OB /9 7 ~ 5, the average solid angle for off-beam 
emission is SIob = 47r[(l — cos6> OB ) + (7r/2 — 8 ob ) sin# OB ] = 
3.4 sr. This gives a limiting distance for EGRET detec- 
tion of off-beam pulsars as du m ~ 140 pc, about halfway 
through the Gould Belt. Assuming the supernova rate in 
the Belt derived by Grenier (2000) of 75-95 Myr^kpc- 2 , 
there would be about 23 - 29 neutron stars of age < 5 
Myr within 140 pc of the Sun, about 8 of which could 
be detected as off-beam pulsars with a beaming factor of 
^ob)/47t ~ 0.27. This estimate is consistent with the 
above relative detectability estimate derived above based 
on geometry. The GLAST out-of-plane sensitivity to point 
sources having a cutoff at 1 GeV is 2 x 10~ 8 phcm~ 2 s _1 
(Digel, private comm.), which would allow GLAST to de- 
tect off-beam pulsars out to a distance of c?i; m ~ 240 pc, 
most of the way through the Belt, and will be able to de- 
tect pulsations of those off-beam pulsars within 150 pc. 
Therefore, GLAST should be able to detect pulsations in 
all of the off-beam pulsars that EGRET detected as point 
sources. 
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Fig. 4. — Model cascade spectra for different observed lines 
of sight, f's. Typical Vela-type pulsar parameters, P = 0.1 s, 
B = 4 X 10 12 G and a = 10° are adopted. The f - a = 4° 
curve is for a typical on-beam pulsar, and the £ — a = 15° and 
f — a = 60° curves are off-beam pulsar spectra. Notice that the off 
beam spectra are intrinsically fainter than the on-beam spectrum, 
and have a high-energy cut off at lower energies. A power law fit 
of the £ — a = 15° spectrum from 100 MeV to 1 GeV gives a spec- 
tral index 2.65 ± 0.57, consistent with the values observed from the 
unidentified "Gould Belt sources". 

3. CONCLUSION 

With the polar cap cascade geometry inferred from the 
known gamma-ray pulsars, we find that there is a large 
phase space for which the line-of-sight misses both the 
main gamma-ray beam and possibly the radio beam, but 
still cuts across a much broader, fainter gamma-ray beam 
produced by the curvature-radiation cooling of the pri- 
mary particles. We identify such off-beam 7-ray pulsars 
as candidates for the new population of the unidentified 
EGRET sources associated with the Gould Belt. The off- 
beam polar cap cascade emission exhibits the low luminos- 
ity and soft spectrum in the EGRET band that is charac- 
teristic of the observed Gould Belt sources. A rough esti- 
mate indicates that the off-beam sources might be ~ 4 — 5 
times more detectable than the on-beam sources due to 
the larger solid angle of the off-beam emission. A more 
accurate prediction of the fraction of off-beam 7-ray pul- 
sars in the Gould Belt will require modeling of the pul- 
sar population and luminosity selection effects. However, 
since EGRET detects all favorably oriented on-beam pul- 
sars and a large fraction of off-beam pulsars in the belt, 
luminosity selection effects may not be a dominant factor. 
GLAST will be able to detect pulsed emission from these 
sources, and the pulse shape should be a broad single- 
peaked profile. 

Off-beam emission is not expected in outer gap mod- 
els (Yadigaroglu & Romani 1995, Cheng & Zhang 1998), 
which predict that most low-latitude (|6| < 5°) uniden- 
tified EGRET sources in the plane are radio-quiet, 
Geminga-like pulsars. Thus, the radio-quiet outer-gap pul- 
sars will be on-average as luminous in 7-rays as the radio- 
loud pulsars. This is in contrast to our prediction that po- 
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lar cap radio-quiet pulsars (at least those detected within 7-rays than radio-loud pulsars. 

a few hundred pc) should be on average less luminous in We thank Isabelle Grenier, Matthew Baring and Scth 

Digcl for valuable input and discussions. 
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